Chaos-induced enhancement of resonant multielectron recombination in highly 

charged ions: Statistical theory 
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A statistical theory of resonant multielectron recombination based on properties of chaotic eigen- 
states is developed. The level density of many-body states increases exponentially with the number 
of excited electrons. When the residual electron-electron interaction exceeds the interval between 
these levels, the eigenstates (called compound states or compound resonances if these states are in 
the continuum) become "chaotic" superpositions of large numbers of Hartree-Fock configurational 
basis states. This situation takes place in some rare-earth atoms and many open-shell multiply 
charged ions excited in the process of electron recombination. Our theory describes resonant multi- 
electron recombination via dielectronic doorway states leading to such compound resonances. The 
result is a radiative capture cross section averaged over a small energy interval containing several 
compound resonances. In many cases individual resonances are not resolved experimentally (since 
the interval between them is small, e.g., < 1 meV, possibly even smaller than their radiative widths), 
therefore, our statistical theory should correctly describe the experimental data. We perform nu- 
merical calculations of the recombination cross sections for tungsten ions W 9+ , q = 18-25. The 
recombination rate for W 20+ measured recently [Phys. Rev. A 83, 012711 (2011)] is 10 3 greater 
than the direct radiative recombination rate at low energies, and our result for W 20+ agrees with 
the measurements. 
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PACS numbers: 34.80. Lx, 31.10.+Z, 34.10.+X, 32.80.Zb 



I. INTRODUCTION 

The majority of atoms in highly excited states and 
some open-shell (e.g., rare-earth) atoms even in vicinity 
of the ground state behave as complex, chaotic many- 
electron systems with very dense spectra and strong con- 
figuration mixing. The many-electron wave function in 
such a system is a mixture of a large number of many- 
excited-electron basis states (Slater determinants) with 
nearly random mixing coefficients. Ordinary theoretical 
methods, such as configuration interaction, which work 
well for atomic systems with few (say, two) electrons 
above closed subshells, become practically useless here. 
On the other hand, such complex systems can be de- 
scribed using statistical approaches. Statistical methods 
for chaotic compound states are widely used in nuclear 
physics (see, e.g., Refs. [H-Hj]). Similar methods were de- 
veloped for open-shell atomic systems [l-Q ■ In Ref. Q we 
used the cerium atom as a testing ground for the appli- 
cations of the statistical theory, and studied properties of 
the Hamiltonian matrix in chaotic many-body systems, 
statistics of the intervals between energy levels, average 
orbital occupation numbers as functions of the excita- 
tion energy, electromagnetic amplitudes between chaotic 
many-body states, enhancement of weak interaction ef- 
fects in such states, and electronic and electromagnetic 
widths of chaotic compound resonances. The statistical 
theory was later also tested for multicharged ions 

In the present paper we want to use the statistical the- 
ory to describe the effect of compound resonances on elec- 
tron recombination with open-shell multicharged ions. 
Electron recombination is an important process in lab- 



oratory and cosmic plasmas, as well as in ion storage 
rings. Theory and experiment agree very well for rela- 
tively simple systems with one or two valence electrons 
above closed shells (see, e.g., Ref. [l(| )■ For more com- 
plex systems theory and experiment often deviate signif- 
icantly. For example, strong enhancement of the recom- 
bination rate were observed for Au 25+ [11], U 28+ [l2[ 
and W 20+ [13[ at low electron energies. For such ions 
the observed recombination rate is orders of magnitude 
greater than that due to direct radiative recombination 
(RR) . In simpler systems the enhancement is due to reso- 
nant dielectronic recombination, though even for ions like 
Fe 9+ 3p 5 and Fe 9+ 3p 4 the dielectronic recombination ap- 
pears to be deficient [l4| . In complex open-shell ions such 
as Au 25+ 4/ 8 and isoelectronic W 20+ the absolute major- 
ity of the resonances correspond to many-excited-electron 
eigenstates which have very high density (exponentially 
small level spacings). In Refs. we used a statistical 
approach to show that the 200-times enhancement over 
RR observed in Au 25+ [ll| is due to electron capture in 
these compound resonances. 

In this paper we use the statistical theory to calculate 
the recombination rates for tungsten ions W 9+ , q = 18 
25. Our results for W 20+ , where the measured rate at low 
(~ 1 eV) energies is 10 3 times higher than the radiative 
rate [l3l ]. are in good agreement with experiment. For 
other tungsten ions we predict similar strong enhance- 
ments of the recombination rate. 

A detailed derivation of the statistical theory for the 
recombination cross section is presented in the next sec- 
tion. Numerical calculations based on our statistical the- 
ory are somewhat similar to those of standard dielec- 
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FIG. 1: Schematic diagram of electron capture into a strongly mixed multiconfigurational eigenstate (dark block) through a 
dielectronic doorway configuration h" 1 ^. Dashed lines show the Coulomb interactions. 



tronic recombination. However, in the statistical theory 
one does not need to diagonalize potentially very large 
Hamiltonian matrices for the excited states of the com- 
pound ion. Instead, the statistical formalism contains 
new parameters, namely the spreading width T spx and 
average electron orbital occupation numbers in the ini- 
tial (ground) state of the ion. The total cross section of 
electron capture into compound resonances is given by 



7T 



r 

abh,lj 



\(a,b\v\h,elj) - (b,a\v\h,elj)\ 2 



{e-e a -e b + e h ) 2 + rf pr /4 
x (n h h c (l - n a )(l - hb))i (1) 



where e and k are the energy and wavenumber of the 
incident electron, (a,b\v\h,elj) and (b, a\v\h, elj) are the 
direct and exchange Coulomb matrix elements, respec- 
tively, the sum is over the single-particle states (orbitals) 
of the hole (h) and excited electrons (a and 6), and the 
partial waves Ij of the incident electron (c), and n a , hj,, 
etc., are the corresponding occupation numbers. The ex- 
pectation value (n^n c (l — fi a )(l — hb))i over the initial 
target state tells one that to transfer electrons from or- 
bitals h and c into a and b, the former must be at least 
partially occupied and the latter at least partially empty. 
The continuum electron wave function c is normalized to 
the delta- function of energy and h c = 1. 

The presence of the spreading width r spr in Eq. ([1]) 
is due to the fact that the two-electron-one-hole exci- 
tation h~ 1 ab is not an eigenstate of the highly excited 
ion. This state is embedded in the dense spectrum of 
other, multiply-excited states and is strongly mixed with 
them. This dielectronic state plays the role of a doorway 
state which "decays" to more complicated states, and the 
width of this "internal" decay is denoted by r spr . This 
process is faster than either autoionization or radiative 
decay ("external" decays), and r spr is greater than the 
autoionization or radiative widths (see below). 

FigureQ]presents a perturbative, diagrammatic picture 
of electron capture through the doorway /i -1 afr. In this 
temporal picture the process looks as a series of electron 
collisions. The initial electron c collides with an atomic 
electron in state h and excites it into state b. Then one of 
the excited electrons interacts with another atomic elec- 
tron to produce more excitations, etc. As a result, the 
initial electron energy is shared between many electrons, 
and none of them has enough energy to escape. In this 



way a long-lived compound resonance is formed. A sim- 
ilar picture of neutron capture by nuclei dates back to 
Niels Bohr [l5j and is very well known to nuclear physi- 
cists. This temporal picture assumes that individual time 
steps St can be resolved. However, the uncertainty rela- 
tion, 8t5E > h, would then require a large energy un- 
certainty 8E. In the recombination process the total en- 
ergy of the system "electron + ion" is well defined. This 
means that all the components (steps in the process in 
Fig.[I| are present in the long-lived quasistationary com- 
pound state which captures the electron. The language 
of strong configuration mixing is more appropriate in this 
case. With a perfect energy resolution one would see a 
very dense spectrum of narrow, possibly overlapping res- 
onances. Broad doorway dielectronic states (with width 
r spr ) can only introduce a variation of the average height 
of these narrow compound resonances on the energy scale 
Ae ~ r spr . 

Compared with the total resonant capture cross sec- 
tion, the recombination cross section 
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contains an additional factor ujf, known as the fluores- 
cence yield. It accounts for the probability of radiative 
stabilization of the resonances (as opposed to autoioniza- 
tion) , and is given by 
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where and r( a ) are the resonance radiative and au- 
toionization widths, respectively. Expressions for these 
widths as well as the spreading width r spr , are presented 
in the next section. 

Note that the capture cross section (fTJ) is not very sen- 
sitive to the specific value of the spreading width r spr 
which for multicharged ions is about 0.5 a.u. [7[ (see Ta- 
ble HI . After angular reduction of the Coulomb matrix 
elements in Eq. (TTJ) [Sec. HH Eq. (|30|) ]. numerical calcu- 
lations of the capture cross section are straightforward. 

An additional simplification occurs in heavy open-shell 
ions like Au 25+ and W 20+ , which have almost unit fluo- 
rescence yield. Indeed, the compound states in these ions 
contain very large numbers of principal basis-state com- 
ponents, N ~ 10 4 . Each component contributes to the 
radiative decay into a large number of states below this 
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compound state. In contrast, only one or few dielectronic 
(doorway) state components have nonzero Coulomb ma- 
trix elements that allow electron autoionization into the 
continuum (see Fig. [T] and Eq. ((IJ). Therefore, the au- 
toionization width of the compound resonance is sup- 
pressed by the small weight factor 1/N of the dielectronic 
components in the compound state. This means that the 
captured low-energy electron cannot escape, i.e., after 
the capture the radiative process happens with nearly 
100% probability. (A similar effect in neutron capture 
by nuclei is described, e.g., in Ref. [Hj]). In this situation 
p( Q ) <§; pO) anc i w ^ ~ p so that the electron recom- 
bination cross section, Eq. ([2]), is independent of the 
radiative width. In this regime one observes maximum 
chaos- induced enhancement of the resonant multielectron 
recombination. For example, the electron capture cross 
section calculated using Eq. (p} for Au 25+ [8J, was found 
to be in good agreement with experiment at low energies. 
On the other hand, in ions with a smaller number of ac- 
tive electrons, the number of components N may not be 
so large, leading to r( a ) > r' r ) and cof < 1. 



II. THEORY 

A. Resonant recombination cross section 

The resonant radiative electron-ion recombination 
cross section is given by the sum over the resonances v 
with the angular momentum and parity J n , 



(4) 



where g(J) = (2J+l)/[2(2Jj+l)] is the probability factor 
due to random orientation of the electron spin and angu- 
lar momentum Ji of the target ion (l8| . The individual 
resonant contributions cr„ are given by the Breit-Wigner 
formula [181 ] . 
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where Tu is the the radiative decay rate (the total "in- 
elastic width" ) , r^f' is the autoionization decay rate (the 
"elastic width"), and T u = + rl a) is the total width 
of the level v. (We assume here that other inelastic chan- 
nels, e.g., electronic excitation, are closed, which is cor- 
rect at low incident electron energies) . The energy of the 
i/th resonance is given relative to the ionization threshold 
/ of the final compound ion, e v = E v — I. 

For systems with dense compound resonance spectra 
the recombination cross section displays rapid energy de- 
pendence, that may not even be resolved experimentally. 
Thus it is natural to average the cross section over an en- 
ergy interval Ae which is large compared with the small 
mean level spacing Dj^ and the total resonance width, 



but much smaller than e. This gives 
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where the integration limits are formally infinite, since 
the contribution of each resonance to cross section de- 
creases rapidly away from e ~ £„. The number of reso- 
nances with a given J 71 ' within Ae is Ae/Dj^, and after 
averaging, the recombination cross section Q becomes 
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Here (...) means averaging of the width factor at the 
given energy. 

If the fluorescence yield w/ = r^ r ^/r fluctuates weakly 
from resonance to resonance, the recombination cross sec- 
tion ay can be factorized, i.e., a r c — uifa c . For uif w 1 
the energy-averaged capture cross section 



(2J + 1) (r(°)) 

k 2 ^ (2 Ji + 1) Dj* 



(8) 



is the same as the recombination cross section. 

In the opposite case of small radiative widths, autoion- 
ization dominates (r (r) < T (q) ) and w/ <C 1, so Eq. (jTj) 
yields the recombination cross section in the form 



- 2 < r(r) ) y (2J 
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where pj* — \/Dj* is the level density and (r( r )) is 
given by Eq. ([3"I|. Note that the sum in Eq. © is 
the total density of states, which can be found without 
constructing states with definite J. 

To estimate the recombination cross section in the gen- 
eral case one can use the following formula, 



(10) 



It follows from Eqs. ([T])-© if w/ does not depend on J. 



B. Nature of chaotic compound states 

The density of excited states p(E) in a many-electron 
ion, especially with an open shell, increases rapidly (ex- 
ponentially) as the number of excited electrons increase. 
Consider n electrons that can be distributed among a 
number of single-electron states g — J^i 2(2Z + 1), where 
I is the orbital angular momentum of the subshells avail- 
able. The total number of many-body states that can be 
constructed is given by 



expfn ln(g/n) + n] 



n\(g — n)\ 



(11) 



4 




(a.u.) 



FIG. 2: Level density in Au 24+ . The black dots are the result 
of the numerical calculation [t| . The solid line is the analytical 
fit, p(E) — AE~ V exp(ay/E), motivated by the level density 
calculated using the Fermi gas model 0, @| . The inset shows 
the densities of states with different J at the ionization energy 
E = I. 



where we used the Stirling formula and assumed g 3> n. 

Equation (1111) indicates exponential increase of the 
number of many-electron states and the corresponding 
decrease of the energy interval between them as the num- 
ber of "active" electrons n increases. For example, Fig. [2] 
shows how the density of multielectron excited states of 
Au 24+ increase with energy E. The small level spac- 
ings between the states mean that even a small resid- 
ual electron-electron interaction will cause strong non- 
perturbative mixing of the many-electron configuration 
basis states (Slater determinants) \$k)- This occurs when 
the off-diagonal matrix elements of the Hamiltonian Hij 
become greater than the energy spacing Dy between the 
basis states i and j coupled by the residual interaction, 



When the mixing is strong, each eigenstate 



£<#°i**>- 



(12) 



contains a large number N of principal components 
These are the basis states for which the expan- 
sion coefficients (which satisfy the normalization condi- 



tion £ fc \C 
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1) have typical values CY 
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The number of principal components in an eigenstate can 
be estimated as N ~ T spr /D, where 



F 



spr 



D 



(13) 



is the spreading width, and D is the mean level spacing 
between the basis states (or eigenstates) . Such eigen- 
states are called compound states, and are well known, 
e.g., in nuclear physics literature [l|. 
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FIG. 3: Components of the 590th eigenstate with J n — 13" /2 
in Au 24+ from a two-configuration calculation [3] and C\ (E) 
(histogram) fit by the Breit-Wigner formula, Eq. (I14j> (solid 
line). The two configurations, 4/5/ 2 4/iy 2 5pi/ 2 5p3/25/ 7 /2 and 
4/5/24/7/25pi/25d 3 / 2 5(;7/ 2 , produce a total of 143360 many- 



electron states with J from ~ 



tof. 



For example, in Au 24+ the mean spacing between the 
excited states with a given angular momentum and par- 
ity, near the ionization threshold is Dj^ ~ 1 meV and 
r spr - 10 eV, so that N - r spr /L> ~ 10 4 0. Numer- 
ical calculations involving a relatively small number of 
configurations confirm that in this case the eigenstates 
are indeed chaotic superpositions of the basis states, see 
Fig. El 

The energies Ek of the principal basis components lie 
within the spreading width of the eigenenergy E v of the 
compound state, \Ek — E v \ < r spr . The components out- 
side the spreading width decrease quickly, so that they 
do not give much contribution to the normalization. It 
was tested in Refs. [1, 0] that components of the chaotic 
eigenstates have the statistics of Gaussian random vari- 
ables with zero mean (Fig. [31 top). On the other hand, 
the variation of their mean-squared value as a function 
of energy (Fig. [31 bottom) can be approximated by a 
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simple Breit-Wigner profile, 



r 2 /4 

x spr/ 



I 2 = AT" 1 

1 k 1 {E k ~E v y+vi v j^ 

with TV = 7iT S p r /2Z? fixed by normalization 



(14) 
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\c[ u) \ 2 dE k /D=l. 



To summarize, the chaotic compound states have the 
following properties: (i) Each eigenstate contains a large 
number N of principal components Ci ~ 1/s/N, cor- 
responding to the basis states \$> k ) which are strongly 
mixed together, (ii) Owing to the strong mixing, the 
only good quantum numbers that can be used to classify 
the eigenstates, are the exactly conserved total angular 
momentum and parity J* and the energy, (iii) The de- 
gree of mixing in this regime is in some sense complete, 
i.e., all basis states that can be mixed (within a certain 
energy range) are mixed. The notion of configurations 
based on the single-particle orbitals becomes largely ir- 
relevant for the purpose of classifying the eigenstates. 
Each eigenstate contains substantial contributions of all 
nearby configurations. 

These properties of chaotic compound states enable 
one to calculate the mean-squared matrix elements of dif- 
ferent operators without diagonalization of prohibitively 
large configuration-interaction Hamiltonian matrices. 



C. Mean-squared matrix elements between 
compound states 

Consider a two-body operator (e.g., the Coulomb in- 
teraction) 



V = ~ y^ y (ab\v\hc)alala h a c 



abch 



A matrix element of V between two compound states, 
and is given by (see Eq. (fH)) ) 

(9 v \V\*i) =^Ci v >C^^ k \V\^), (15) 



or 



(*»\V\%) = 1^2 i(*b\i\hc) - (ba\v\hc)) 

(ab)ch 

x {9 v \alala h a c \^i), 



(16) 



where in the last equation the sum is restricted to dis- 
tinct pairs {ab), and v \a\a\aha c \^> \) determines the 
contribution of the two-particle transition ch —¥ ab. Due 
to the assumption that the expansion coefficients for 
chaotic compound states are random and uncorrelated 



{C { k v) = C^C^ =0iorv^i), the value of the matrix 
element averaged over many compound states v is zero, 



(* I/ |V'|* i ) = 0, {^ v \a\a\a h a c \^ i )=§, (17) 



To determine the autoionization width (Sec. Ill D[) we 
need to calculate the mean-squared matrix element. It is 
derived using the statistical properties of the expansion 

coefficients, C { k v) *c\ v) = \C { k u) \ 2 6 kl , so that 



{ty v \aUla c ah\Vi)(Vi\alal,a b ,a a ,\V 



t J 



o'aa>o'bb'0 : cc'0'hh>\( x S'v\alala h a c \ ll S'i)\ 2 . (18) 
Hence, the mean-squared matrix element is 

l(*,im>l 2 = \"£\(ab\e\kc)-{ba\v\hc}\ 3 

abch 



x \{y v \alala h a c \*i)\*. (19) 
Let us introduce the strength function 



w{E k ;E v ,T spI ,N)=C { k )2 , 



(20) 



which describes the spreading of the component k over 
the eigenstates v {C^ are assumed to be real). This 
function depends on the number of principal compo- 
nents N of the eigenstate, Eq. (|12[) . the spreading width 
r spr , and the difference E v — E k between the energies of 
the compound state and component k. In the simplest 
model (1J w(E k ; E u ,T spr , N) is a Breit-Wigner function, 
cf. Eq. ([TU). Using Eq. ([T5|). we then obtain 



\{MaUWa c \^i)?=Y,Y. C ^ C i C * C v 

kk 1 IV 

x {§ k \a\ala h a c \® k '){®l< \ala{a b a a \^i) 

kk' 

x{^ k ,\ala\a b a a \^ k ){^ k \ala\a h a c \^ k ,). (21) 

To obtain the last expression we used the properties 
of the components and the definition (I2TJ1) . and de- 
noted Wi(E k >) = w(E k r,Ei,ri%,Ni) and w u (E k ) = 

w(E k ;E v ,rj&,N v ). 

We can assume, without the loss of generality, that 
the number of principal components \3> k ) in state v is 
greater than or equal to the number of components \<& k ') 
of state i, i.e., ripl/D^ > ri p \/Di. The matrix ele- 
ment (& k \a'l l ala} l a c \Q k ') does not vanish only if \$ k ) = 
a\,a\a h a c \® k i) , so that E k -Ey ~ e a +e b -e h -e = u) abtCh . 
Using closure to sum over k in Eq. (|21[) . we obtain 



K^^alalahad^i)^ = ^Wi(E k >)w v (E k > + uj a b,ch) 

k' 

x{$ k ,\n h n c (l - n Q )(l - h b )\<S> k ,).(22) 



6 



In deriving this equation we used the anticommutation 
relations satisfied by the creation and annihilation op- 
erators, and introduced the the occupation number op- 
erators n a = a^ a a a . The matrix element (<&w\h}Ji c (l — 
h a )(l — h b )\<$>k') is equal to unity if the orbitals h and c 
are occupied, while the orbitals a and b are vacant in the 
state i.e., the transition ch — > ab is possible. 

If one assumes that the single-electron-state occupan- 
cies vary slowly with the excitation energy, then the ma- 
trix element of the operator ra^n c (l — n a )(l — h b ) in 
Eq. ([22]) can be replaced by its expectation value, 

^w t {E k ,)(<$> k ,\h h h c (l - h a )(l - h b )\$k>) = 

k' 

(h h n c (l-h a )(l-n b ))i, (23) 

subject to the normalization condition Y)u Wi{E k i) = 1. 
The right-hand side of Eq. (|23[) is the value of the oc- 
cupancy times "emptiness" in the compound state 
averaged over a number of neighboring states. 

Replacing the matrix element (3?fc'| . . . \*&k') by its av- 
erage (l23l) in Eq. (|22[) . and changing summation to inte- 
gration, one obtains 



\{^ uWaa\a c ah\^ i)\ 2 = (n h h c (l - n a )(l - h b )) l 



X / Wi{E V )w v {E k > +UJ a b,ch) 



dE k , 



This result can be written in the following form: 



(24) 



\{^ u\a\a\a h a c \^ i)\ 2 = (n h n c (l - n a )(l - n b )) l 



xW(rW,rM,A). (25) 



In this expression 



pr' spr? / / w V « j -^t? 1 - spr 

dE 

xw(E kl +LJ ab , ch ;E u ,T^ t ,N u )—^, (26) 



where A = E v — Ei — uj ab ,ch, is a "spread" delta- function 
which was studied in Refs. 0,110. It peaks at A = 
and describes the approximate energy conservation for 
the transition between compound states induced by the 
two-electron transition ch —¥ ab and broadened by the 
spreading widths. For the Breit-Wigner strength func- 
tions one has 



S( T (i) r M A ) = 1 Tspr 

spr' spr' '—^i 



2tt A^+r^' 



(27) 



where T spr = ri% + T { S pl From Eqs. (H]) and ((25j>, the 
mean-squared matrix element of the two-body operator 
between the compound states is finally obtained as 

IWW=J E \{ab\v\hc)-{ba\v\hc)\ 2 

(ab)ch 

x(n h h c (l-n a )(l - n b ))iDj(TW,r%l,&). (28) 



In this expression the summation is carried out over the 
single-electron states a, b, h, and c. Note that if l^j) 
is a simple, unmixed state, there is no sum over k' in 
Eqs. In this case T { s % = and r spr = T%\ in 



Eq. d27J). 

For a one-body operator M = J2 a b M ab a\a b Ahe mean- 
squared matrix element is obtained similarly 0, H H| : 



\(* v \M\* f )\2 = ^\{a\m\b)\\n a (l-n b )) v 



ab 



xDfSF&r&Ef-Ev-Uba), (29) 

where uj ba = e b — s a is the energy of the single-electron 
transition a — > b. 



D. Capture cross section and autoionization width 

The autoionization width rl a) = 2ir\(^ „\V\^i)\ 2 gives 
the transition rate between the initial state, e~ + A q+ , 
and the multiply exited compound resonance of the ion 
^4(9-1)+ to the two-body Coulomb interaction V. 
Unlike the complex multiply excited states | X I , ; ,), the ini- 
tial state of the recombination process is simple. It de- 
scribes an electron with the energy e incident on the 
ground (or low-lying excited) state |$o) 01 the target, 
which is often dominated by one configuration. It is 
clear that the autoionization width averaged over com- 
pound resonances is determined by the mean-squared 
matrix element of the Coulomb interaction between elec- 
trons, which is given by Eq. (f2"5|). The initial state 
l^i) = 1^0) c ) is thus a compound state with negligi- 
ble spreading width Tspr <C ripl. The total width of 
the function 5, Eq. (|27j) . is dominated by the compound 



resonance width T 



spr 



r 



O) 

spr ■ 



Nonzero contributions to the sum in Eq. (|28|) . i.e., to 
V( a \ arise from the basis states which differ from the ini- 
tial state |$o, c) by the single-particle states of two elec- 
trons. Therefore, it is sufficient to sum over the doubly- 
excited basis states, 



1*0=2*5^ 



M2 



\(<S> d \V\<S> ,c)f 



Such two-electron excitations |<&d) play the role of door- 
way states for the electron capture process. Since these 
states are not the eigenstates of the system they have a 
finite energy width r spr . The wave function of a doorway 
state can be constructed using the creation-annihilation 
operators, = a\ a\a h \^o, c), where a = n a l a j a m a 

and b = n b l b j b m b are excited single-electron states, and 
h = nhlhjhmh corresponds to the hole in the target 
ground state. Of course, to form the doorway states with 
a given total angular momentum J, the excited electrons 
and the ionic residue must be coupled into J. However, 
the 2 J + 1 factor and summation over J in Eq. ([5]) ac- 
count for all possible couplings. This means that the sum 
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over the eigenstates in Eq. ([8j can be replaced by the sum 
over the one-hole-two-electron excitation, as in Eq. (|28|l , 
and one obtains the capture cross section in the form of 
Eq. p. 

Note that when the number of active electrons and 
orbitals is large, the occupation numbers for different 
orbitals become statistically independent. In this case, 
the correlated product of the single-particle occupancies, 
Eq. (j2"3")l can be approximated by the fractional occupa- 
tion numbers of the electronic subshells with definite j. 
The orbital c is taken a continuum, c = eljm in Eq. (JlJ. 
Its wave function is normalized to the delta-function of 
energy, and it is occupied in the initial state, i.e., h c = 1. 
After summation over the magnetic quantum numbers 
m a , mb, etc., and angular reduction of the Coulomb ma- 
trix elements, the final expression for the capture cross 
section is 



abhdj y 



spr 



E 



e a -e b + e h ) 2 + T 2 pi /A 
a,b\\V\\\h,elj) r 



2A + 1 

A 

£(-l)A+V+l 



(a,b\\V x \\h,slj)-(2X+l) 



A ja 3 
A' jb jh 



(b,a\\Vy\\h,elj) 



2jh + 1 



2ja + 1 



rib 



2jb + 1 



(30) 



Here n a , rib and rih are the occupation numbers of the 
corresponding subshells (ranging from to 2j a + 1, etc.), 
and e a , £b, and Sh are their energies. The two terms in 
square brackets represent the direct and exchange con- 
tributions, and (a, b\\V\\\h,elj) is the reduced Coulomb 
matrix element: 



(a,b\\Vx\\h,c) = ^(2j a + l)(2j b + l)(2j h + l)(2j c + 1) 
x Z(l a + l c + \)Z(l b + l h + \) (31) 



A ja jc\ /A j b j h 



0-H/lO 

. u 2 2/ V 22 

where £(L) = [1 + ( — l) L ]/2 is the parity factor, and 



1 j Rx(a,b;h,c), 



R\(a,b;h,c) = 



r [/-W/c(r)+5a(r) flc (r)] 



x [f b (r')f h (r') + g b (r r )g h (r')}drdr' (32) 

is the radial Coulomb integral, / and g being the upper 
and lower components of the relativistic orbital spinors. 

Once a c is known, Eq. (|8j allows one to estimate the 
average ratio T^/D for a typical J*, 



r(°) 



k 2 (2Ji + l)a c fc 2 (2J, + l)a c 



(33) 



where the sum in the denominator is over the angular 
momentum and parity J 7 ' which contribute effectively to 



the capture cross section. For example, J max ~ 10 and 
J. i = 6 for the recombination of Au 25+ and W 20+ . A 
typical distribution of level densities pj^ for different J, 
is shown on the inset of Fig. [2J 



E. Radiative width 

The second step of the recombination process is radia- 
tive stabilization. Any excited electron in the compound 
state |\I/„) can emit a photon. Using Eq. (|21)|) the total 
photoemission rate T^ r ' can be estimated as a weighted 
sum of the single-particle rates, 



4w 



ba 



^ 3c 3 

a,b 



\(a\\d\\b)f 



2j b + 1 



1 - 



2j a + 1 



(34) 

where U)b a = — e a > 0, (a||d||6) is the reduced dipole 
operator between the orbitals a and 6, and (. . . ) v is the 
mean occupation number factor. The mean subshell oc- 
cupation numbers for a given energy can be obtained by 
averaging over the basis states involved, e.g., 



n a {E) 



(35) 



(k) 

where n a is the occupation number of the subshell a in 
the basis state k. 

Since have large numbers of principal components 
N, the fluctuations of their radiative widths are small, 
~ 1/ %/iV. This can also be seen if one recalls that a 
chaotic multiply excited state is coupled by photoemis- 
sion to many lower-lying states, and the total radiative 
width is the sum of a large number of (strongly fluctuat- 
ing) partial widths. A similar effect is known in com- 
pound nucleus resonances in low-energy neutron scat- 
tering p], [lj|. In multicharged ions with dense spec- 
tra of chaotic multiply-excited states, the autoioniza- 
tion widths are suppressed as r( a ) oc 1/iV. Physically 
this happens because the coupling strength of the two- 
electron doorways state to the continuum is shared be- 
tween many complex multiply-excited eigenstates. The 
radiative width does not have this suppression since all 
components of a compound state contribute to the radia- 
tive decay. As a result, the radiative width may dominate 
in the total width of the resonances, 

r« > r( Q ), mak- 
ing their fluorescence yield close to unity. Our numerical 
results for the recombination of Au 25+ presented in Q, 
supported this picture. 



III. NUMERICAL RESULTS 



In this section we apply our theory to calculate the 



17^ 



to 



recombination rate for the tungsten ions from W 
W 24+ . Experimental data are available for the recom- 
bination of W 20+ forming W 19+ We will use this 
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1 2 3 

Orbital angular momentum 

FIG. 4: Energies of the occupied and vacant single-particle 
orbitals of W 19+ obtained in the Dirac-Fock calculation. 



with energies close to the ionization threshold, which act 
as the doorway states. One also needs to estimate the 
spreading width r spr . For low-energy electron recombi- 
nation, we restrict the consideration to the energy inter- 
val, 



E = I±AE/2, 



(36) 



where E = for the ground state of the final-state ion, 
and we choose AE ~ r spr . In practice, we start from 
some initial estimate of the spreading width, and sub- 
sequently find a more accurate value using an iterative 
procedure. 

The spreading width is found from Eq. (|13[) by con- 
structing the Hamiltonian matrix for a number of con- 
figurations within the interval (|36[) . and averaging the 
squared offdiagonal matrix elements of the Hamiltonian 
matrix H; 



system as an example to describe the calculations. Cal- 
culations for other ions are similar. 

When an electron recombines with W 20+ , it can be 
captured into an excited state of the compound W 19+ 
ion. Its ground state belongs to the Is 2 . . . 4/ 9 config- 
uration. Figure 0] shows the energies of its relativistic 
orbitals nlj obtained in the Dirac-Fock calculation. All 
orbitals below the Fermi level, Is to 4/, were obtained 
in the self-consistent calculation of the W 19+ ground 
state. Each of the excited-state orbitals above the Fermi 
level: 5s, 5p, etc., was calculated by placing one electron 
into it, in the field of the frozen W 20+ Is 2 . . . 4/ 8 core. 
The energy of the highest orbital occupied (partially) 
in the ground state is £4/ 7/2 = —18.41 a.u. This value 
gives an estimate of the ionization potential of W 19+ : 
I « |£4/ 7/2 | = 18.41 a.u. This value is in agreement with 
NIST data, I = 18.47 a.u. [HI . 

Excited states of the ion are generated by transferring 
one, two, three, etc. electrons from the ground-state or- 
bitals into the empty orbitals above the Fermi level (Fig. 
H|), or into the partially occupied 4/ orbitals. We are 
interested in the excitation spectrum of W 19+ near its 
ionization threshold. This energy (^20 a.u.) is sufficient 
to push up a few of the nine 4/ electrons, and even excite 
one or two electrons from the Ad orbital. However, the 
preceding 4p orbital is already deep enough to be con- 
sidered inactive. Thus, we treat W 19+ as a system of 19 
electrons above the frozen Kr-like Is 2 . . . 4p e core. Note 
that in constructing the excited state configurations, we 
disregard infinite Rydberg series which correspond to the 
excitation of one electron in the field of W 20+ . Rydberg 
states belong to a single-particle aspect of the e~ + W 20+ 
problem, and are not expected to contribute much to the 
recombination cross section in this system. 

Assuming that the fluorescence yield is close to unity 
(see below), we calculate the recombination cross section 
from Eq. (|30l) . Before using this formula, one needs to 
obtain a list of two-electron-one-hole excitations of W 19+ 



H 2 



N S (N S - 1) 



i<j 



(37) 



where the average is taken over N s basis states whose 
energies Ek = Hkk lie within the energy interval (|3"6"j) . 
The level spacing D in Eq. p^|) is found as average energy 
interval between the states 



D = AE/N S 



(38) 



The list of two-electron-one-hole excitations /i _1 a6 which 
contribute to the sum ([30]) is found by checking which of 
these configurations contribute basis states to the sum in 
Eq. (|37f . It is known that the spreading width is a robust 
characteristic of the system. Indeed, we have checked 
that when more configurations are included, both D and 
Hfj decrease, whereas r spr does not change much (see 
also Ref. @). 

When finding Hfj and D we use basis states with def- 
inite projection of the total angular momentum J z cor- 
responding to the minimal value of J z (0 or 1/2), rather 
than the states with definite total angular momentum 
J. This method is significantly simpler than the use of 
the basis states with definite J and J z , and produces the 
same results for r spr , Eq. (fl"3)) . 

Table Q] shows the spreading widths for the compound 
ions of tungsten, W^ 9 " 1 ^, with excitation energies close 
to the ionization threshold, formed in the process of low- 
energy electron recombination with W^ + . With the ex- 
ception of the target ion with the smallest number of 
4/ electrons (W 25+ 4/ 3 ), the spreading widths are in the 
range 0.5-0.7 a.u. In fact, the value of r spr does not 
strongly affect the magnitude of the capture cross section, 
Eq. (|3H)) . since the area under the Breit-Wigner contour 
corresponding to each doorway h~ l ab is independent of 
r 

spr • 

As discussed in Sec. Ill Bt strong mixing of the basis 
states results in the eigenstates with large numbers of 
principal components, N ~ r spr /Z? ~ Hf-/D 2 1. This 
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TABLE I: Electron capture cross sections a c and rate co- 
efficients a c for the tungsten ions W M+ with the open 4/ 
subshell, and properties of the compound ions at 
excitation energies close to the ionization threshold I. 



Target 



r 



D 

10" 4 a 



k r 



pr 

a.u. 



10 



-16 , 



10~ 7 cm 3 /s 



W 18+ 4/ 10 

W 19 +4/ 9 

W 20+ 4/ 8 

W 21+ 4/ 7 

W 22+ 4/ 6 

W 23+ 4/ 5 

W 24 +4/ 4 

W 25+ 4/ 3 



15.5 
17.0 
18.5 
20.0 
21.8 
23.5 
25.2 
27.0 



0.2 
0.1 
0.1 
0.1 
0.2 
0.4 
1.3 
11 



70 
93 
105 
96 
76 
48 
25 
5 



0.56 
0.65 
0.68 
0.68 
0.65 
0.59 
0.50 
0.16 



25 
29 
30 
34 
16 
11 
19 
12 



1.5 
1.7 
1.8 
2.0 
0.98 
0.67 
1.1 
0.7 



"Ionization energy of the final-state ions, Ref. 11911 . 
'Capture cross section from Eq. (130 P and rate coefficient for inci- 
dent electron energy e = 1 eV. 




occurs when 



0.1 1 
Electron energy (eV) 

FIG. 5: Recombination rate coefficient of W 20+ . Dashes line 
is the direct radiative recombination rate, Eq. (|40|) ; solid 
line is the capture rate calculated using the present theory, 
Eq. (|30[) ; dot-dashed line shows the same, taking into account 
the velocity distribution of the electron beam (see text); solid 
circles is the measured rate coefficient 11311. 



K = \]Hl/D^>l. (39) 

Table fl] shows that this criterion is fulfilled for all the 
ions studied, and that the expected number of principal 
components is indeed large, N ~ 10 4 . As explained in 
Sec. Ill E| in this case one can expect large fluorescence 
yields, w/wl, This means that the recombination cross 
section will be at the limit given by the total electron 
capture cross section, Eq. (|30j) . 

In the present calculations of a c , Eq. (|30l) . we also in- 
clude in a semiempirical way the effect of screening of 
the Coulomb interaction between valence electrons by 
core electrons. This is done by introducing the screening 
factors f\ in the two-electron Coulomb integrals, assum- 
ing that these factors depend on the Coulomb integral 
multipolarity A only. The factors were calculated to be 
fi = 0.7, f 2 = 0.8, f 3 = 0.9 0. Coulomb integrals of 
other multipolarities are not modified. The above values 
of the screening factors were found in the calculations for 
other atomic systems. However, in practice they change 
little from one atom to another. 

To compare with experiment for W 20+ [l3j], the cross 
section obtained from Eq. (|30| is converted into the rate 
coefficient a c — a c v, where v is the velocity of the in- 
cident electron. The result is shown in Fig. [5] by the 
solid line. Since the sum in Eq. (|30l) has a weak depen- 
dence on the electron energy, the capture cross section at 
low energies is proportional to 1/e, and the correspond- 
ing rate coefficient behaves as a c oc 1/v. The calculated 
rate agrees well with the experimental data in the energy 
range of 0.1-1 eV. At higher energies the experimental 
rate coefficient tends to drop faster than 1/v. 

Figure [5] also shows the rate coefficient for the direct ra- 
diative recombination. The latter is estimated using the 
Kramers formula for the radiative recombination cross 



section [171 ] (in atomic units), 

^^Sm(4:V (40) 

d 3\/3c 3 fc 2 \n Q kJ K ' 

where Zi is the ionic charge Zi (e.g., Zi — 20 for e~ + 
W 2a+ ), and no is the principal quantum number of the 
lowest unoccupied ionic orbital (n = 5 for W 20+ ) 0). 
The energy dependence of this cross section is close to 
1/e, and the corresponding rate coefficient (dashed line 
in Fig. [5J is three orders of magnitude smaller than the 
measurement in the energy range shown. 

Below e = 0.1 eV the measured recombination rate 
coefficient can be affected by the velocity distribution of 
the electron beam, which is characterized by two temper- 
atures, T|| = 0.15 meV and T± = 10 meV pjjj]. Taking 
this into account (see Eq. (18) in Ref. [8[) reduces the cal- 
culated resonant capture rate below 50 meV (dot-dashed 
line in Fig. [5j , bringing it into closer agreement with ex- 
periment. 

As discussed above, the capture cross section has a 
simple 1/e energy dependence at low electron energies. 
Hence, in Table U we show the cross sections and rate 
coefficients for W q+ (q = 18-25) calculated at one low 
electron energy, e — 1 eV. We see that the largest cross 
section is predicted for the ion with the half-filled 4/ 
subshell. On the other hand, all the cross sections are 
within a factor of three of each other, and much larger 
than what one would expect from the direct RR process, 
Eq. (|40l). 

Of course, one must keep in mind that compared with 
the capture cross section, the recombination cross section 
contains an additional factor a;/. The fluorescence yield 
may be significantly smaller than u>t = 1 for ions, in 
which the degree of mixing is not as large as it is in the 
compound W 19+ ion. In particular, this may be the case 
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for the ions in which the mixing strength K (see Table 
U) and the number of principle components N is not too 
large. In this case one should regard a c as the upper limit, 
and use Eqs. © and (| 10[) to estimate the recombination 
cross section. 

IV. CONCLUSIONS 

A detailed derivation of the statistical theory of res- 
onant electron capture by many-electron ions has been 
presented. Numerical calculations have been performed 
for a number of tungsten ions with a partially filled 4/ 
subshcll. The calculated rate coefficient for W 20+ is in 
agreement with the measurements at low electron en- 



ergy. The present approach can be used to investigate 
other processes mediated by chaotic, multielectronic ex- 
cited states. 

After completing the present work, we became aware 
of Ref. [21[ which considers dielectronic recombination of 
W 20+ . 
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